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Detailed  input  instructions  are  provided.  A sample  problem  is  solved  to 
illustrate  usage  of  the  program  and  also  to  present  the  output.  The 
general  output  scheme  is  explained  and  the  output  variables  are  defined. 
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ABSTRACT 


A description  is  given  of  a computer  program  which  calculates 
the  incompressible  boundary  layer  flow  and  pressure  distribution 
over  an  ax isymmetric  body.  A brief  outline  is  given  of  the  computa- 
tion method.  Detailed  input  instructions  are  provided.  A sample 
problem  is  solved  to  illustrate  usage  of  the  program  and  also  to 
present  the  output.  The  general  output  scheme  is  explained  and  the 
output  variables  are  defined. 

ADMIN ISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Naval  Sea  Systems  Command  under  the 
Direct  Laboratory  Funding  Program  on  the  Hydrodynamics  of  Very  High-Speed 
Submarines,  Element  Number  62545N.  The  work  was  performed  under  internal 
Work  Units  IS 20 -004  and  1500-200. 

INTRODUCTION 

A computer  program  is  documented  which  calculates  the  incompressible 
boundary  layer  flow  and  pressure  distribution  over  an  axisymmetric  body  in 
uniform  flow  at  zero  angle  of  attack. 

The  computation  method,  which  is  more  fully  described  in  Reference  1, 
is  first  briefly  outlined  here.  Instructions  are  then  provided  on  the 
program,  including  a listing  of  the  input  READ  statements,  the  definition 
of  the  input  variables,  a number  of  comments  on  usage  of  the  program,  and 
the  memory  and  computer  time  requirements  of  the  program.  A sample  problem 
is  presented  to  illustrate  usage  of  the  program.  The  output  for  this 
problem  is  presented  and  the  output  variables  are  defined. 

DESCRIPTION  OF  COMPUTATION  METHOD 

The  program  contains  three  major  calculations.  The  Douglas  Neumann 
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method^’  is  first  used  to  calculate  the  potential  tlow  pressure  distribution 

^Huang,  T.T.,  il.T.  Wang,  N.  Santelli,  and  N.C.  (.roves,  "Propeller/Hull 
Interaction  on  Axisymmetric  Bodies:  lheory  and  Experiment,"  DTN5RDC 

..Report  76-0115  (in  review). 

“Smith,  A.M.0,  and  J.  Pierce,  "Exact  Solution  of  the  Neumann  Problem. 
Calculation  of  Non-Circulatory  Plane  and  Axially  Symmetric  Flows  About 
or  Within  Arbitrary  Boundaries,"  Douglas  Aircraft  Company  Report  IS  26988 
-(Apr  1958) . 

‘Hess,  J.L.  and  A.M.0.  Smith,  "Calculation  of  Potential  Flow  About  Arbitrari- 
Bodies,"  from  I’rogress  in  Aeronautical  Sciences,  \ol . 8,  Pergamon  Press, 
Oxford  and  New  iork  (1966). 


over  the  original  body.  In  this  method,  the  body  is  divided  into  a number 

of  frustrums  of  cones,  each  of  which  has  a constant  source  density  on  its 

surface.  A set  of  linear  algebraic  equations  is  used  to  solve  for  the 

strengths  of  these  source  densities.  Once  the  source  densities  are  known, 

the  pressure  coefficients  everywhere  in  the  flow  field,  including  the  body 

surface,  may  be  readily  determined.  This  initial  pressure  distribution, 

with  velocity  ratio  linearly  extrapolated  in  the  stern  region,  is  then  used 
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by  the  Douglas  Cebeci-Smith  (CS)  differential  boundary- layer  method  ' ’ to 
calculate  the  viscous  flow  over  the  body.  This  method  makes  the  standard 
thin  boundary  layer  assumption  of  a constant  pressure  across  the  thickness 
of  the  layer.  The  integral  wake  relations  given  by  Granville  are  then  used 
to  calculate  the  flow  in  the  wake.  These  relations  are  essentially  based  on 
experimentally  measured  wake  data  behind  a body  of  revolution. 

The  calculated  displacement  thicknesses  from  the  boundary  layer  and 

wake  methods  are  then  used  to  generate  a new  overall  body-wake  displacement 

model.  In  the  stern/near- wake  region,  where  neither  of  the  two  methods 

properly  model  the  thick  boundary  layer,  a fifth-degree  polynomial  is  used 

to  fair  the  calculated  stern  displacement  thickness  surface  into  the 

calculated  wake  displacement  thickness  surface.  The  polynomial  coefficients 

are  selected  to  provide  continuous  displacements,  slopes,  and  curvatures  at 

the  intersections  of  the  polynomial  with  the  calculated  stem  and  wake 

displacement  thickness  surfaces.  The  Douglas  Neumann  method  is  again  used 

to  calculate  the  pressure  distribution  over  the  resulting  body-wake  displace- 

ment  body,  which  is  then  used  to  calculate  the  new  boundary  layer  flow  over 

^Cebeci,  T.  and  A.M.O.  Smith,  Analysis  of  Turbulent  Boundary  Layers, 
rAcademic  Press,  New  York  (1974). 

cebeci,  T.,  G.J.  Mosinskis,  and  A.M.O.  Smith,  "Calculation  of  Viscous  Drag 
and  Turbulent  Boundary-Layer  Separation  on  Two-Dimensional  and  Axisymmetric 
Bodies  in  Incompressible  Mows,"  Douglas  Aircraft  Company  Report  No. 
MPC-J0973-01  (Nov  1970)  . 

cebeci,  T.,  G.  Mosinskis,  and  I,.C.  Wang,  "A  1 inite-Di fference  Method  for 
Calculating  Compressible  Laminar  and  Turbulent  Boundary  Layers,  Part  II 
..User's  Manual,"  Douglas  Aircraft  Company  Report  DAC-67131  (May  19b9)  . 
Granville,  P.S.,  "The  Calculation  of  the  Viscous  Drag  of  Bodies  of  Revolution, 
DTMR  Report  849  (.Jul  1953)  . 


the  body  and  in  the  wake.  This  process  is  repeated  until  successive  pressure 
distributions  agree  to  within  a specified  error  criterion  at  all  points  along 
the  body-wake  displacement  surface  or  until  the  specified  maximum  number  of 
iterations  is  reached.  A sketch  of  the  original  axisymmetric  body,  the  body- 
wake  displacement  surface,  and  the  definition  of  coordinate  systems  is  given 
in  Figure  1 . 

In  order  to  simplify  the  overall  program  and  also  to  reduce  the  number 
of  input  variables,  a number  of  options  in  the  two  Douglas  programs  which  are 
not  directly  applicable  to  the  present  problem  have  been  deleted.  The  major 
deletions  include  the  effect  of  compressibility  contained  in  both  programs, 
the  calculation  of  vorticity  and  cross  flows  in  the  Dougl as-Neumann  program, 
and  the  option  for  calculating  two-dimensional  flows  in  the  Douglas  CS 
program. 

On  the  other  hand,  the  output  of  the  Douglas  CS  program  has  been 
expanded  in  two  areas.  First,  whereas  the  original  program  printed  out  only 
the  tangential  velocity  profile,  the  present  program  also  prints  out  the 
normal,  axial,  and  radial  velocity  profiles  in  the  boundary  layer.  It  should 
be  noted  that  wake  distributions  in  the  propeller  disk  plane  are  usually  given 
in  terms  of  the  latter  two  profiles.  Secondly,  three  different  methods  are 
used  to  compute  the  overall  drag  acting  on  the  body.  Two  of  these  methods, 
by  Squire-Young'’  and  Granville,  are  essentially  empirical  formulas  for  the 
drag  based  on  boundary  layer  parameters  at  the  stern  of  the  body.  They  arc 
commonly  used  in  cases  where  the  actual  pressure  distribution  over  the  body 
is  unknown.  The  third,  and  more  accurate,  method  consists  of  summing  the 
components  of  drag  due  to  friction  and  pressure.  The  resulting  drag 
coefficients  are  given  in  terms  of  three  different  reference  areas:  frontal 

area,  wetted  area,  and  (volume)"^'’. 


INPUT  INSTRUCTIONS 

INPUT  STATEMENTS 

The  input  statements  by  means  of  which  data  are  entered  into  the 
program  are  as  follows: 

READ (5 ,1012)  IPROP,  NXMM,  IPF,  IVF,  FLPRNT,  TITLP 
*REAP(60)  (TX1(I),  TY1(I),  I = 1,  NXMM) 

**READ(64)  (TX1(I),  TY1(I),  I = 1,  NXMM) 

RF.  AD  ( 5 . S 0 1 0 ) NXT,  LGlb , LG17 
READ (5 , 5051)  ROMAX,  RL,  Ul,  RI 
•READ (5 ,1013)  NTOT,  ICUT,  ICP 
‘READ (5,5025)  XA,  XTL,  XWAK,  XSLP,  CPFRR 
*READ(5,2)  (XX(I),  I = NN  + 1,  NTOT) 

**READ( 22)  (UE ( I) , I = 1,  NXMM) 

**READ(5 ,1013)  N 

**READ(5 ,2)  (X(I),  I = 1,  N) 

**READ(5, 2)  (Y ( I ) , I = 1 , N) 

*Skip  when  velocity  distribution  is  input,  i.e.,  IPROP  > 1. 

**Skip  when  velocity  distribution  is  not  input,  i.e.,  IPROP  < 0. 

The  corresponding  FORMA  1 statements  are  as  follows: 

1012  FORMAT  (415,  F8.4,  15A4) 

5010  FORMAT  (14,  711) 

5051  FORMAT  (3F10.4,  F12.2) 

1013  FORMAT  (2014) 

5025  FORMAT  (SF10.4) 

2 FORMAT  (8F10.7) 

DEFINITION  OF  INPUT  VARIABLES 

IPROP  >1  velocity  (pressure)  distribution  is  input 
<0  velocity  distribution  is  not  input 

NXMM  Number  of  body  points  (typically,  NXMM  = 140  to  loO) 

IPF  >1  output  from  potential  flow  program  is  printed  after 

each  iteration 
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IVF 


FLPRNT 

TITLP 

Txicn 

TY1(I) 

NXT 

LG  lb 

LGl? 

ROMAX 
RL 
U1 
R I 

NTOT 

icirr 

icr 


<0  output  from  potential  flow  program  is  printed  only  after 
initial  calculation  on  original  body 

>1  full  output  from  boundary  layer  program  is  printed  after 
each  iteration 

<0  full  output  from  boundary  layer  program  is  printed  only 
after  final  iteration 

Minimum  value  of  X/RL  for  which  velocity  profiles  are  printed 
after  final  iteration 

Title 

Axial  distance  measured  from  nose  in  feet 
Radius  of  body  at  X = TX 1(1)  in  feet 

Index  of  body  station  where  flow  becomes  turbulent;  if  transition 
is  to  be  calculated  by  program, 

NXT  > NXMM 

= 0 if  transition  point  is  not  to  be  calculated  by  the  program 

= 1 if  transition  point  is  to  be  calculated  by  the  program 

= 0 if  transition  is  instantaneous 
= 1 if  transition  is  gradual 

Maximum  radius  of  body/RL 

Reference  length  = body  length  = TXl(NXMM) 

Free-stream  velocity  in  feet/second 

Reynolds  number/ RL 

Number  of  body  points  + number  of  wake  points,  must  be  <_200 

Maximum  number  of  iterations  for  pressure  on  body-wake  displace- 
ment model,  not  including  initial  potential  flow  calculation  on 
the  original  body  (typically  set  equal  to  3) 

= 0 if  pressure  coefficient  from  previous  iteration  is  used  to 
calculate  boundary  layer  flow 

'•l  if  average  of  pressure  coefficient  from  previous  two 
iterations  are  used  to  calculate  boundary  layer  flow 


XA 


Value  of  X/RL  where  initial  guess  for  wake  velocity  reaches 

0.99  of  free  stream  velocity  UI  (typically  set  equal  to 
1.15) 


XTL  Value  of  X/RL  where  initial  potential  flow  velocity  distribution 

is  linearly  extrapolated  to  tail  (usually  set  equal  to  0.95) 

XWAK  Upstream  value  of  X/RL  where  5th-degree  polynomial  matches  with 

computed  displacement  body  (usually  set  XWAK  equal  to  0.95  or 
two  body  stations  ahead  of  the  separation  point,  whichever 
occurs  first) 

XSLP  Downstream  value  of  X/RL  at  matching  point  (usually  set  equal 

to  1.05) 

CPERR  Maximum  allowable  error  in  pressure  coefficient  CP  along  the 

entire  body  and  wake  (typically  set  between  0.005  and  0.02) 

XX(I)  Values  of  X/RL  for  the  wake  displacement  model  (typically 

1 ,0<  XX (I)  < 30) 

UE(I)  Input  velocity  ratio  at  the  input  body  points  TX1(I),  TV  1(1) 

N Number  of  points  where  changes  to  input  velocity  ratio  are  made 

X ( I ) Values  of  X/RL  where  velocity  ratio  changes  are  made,  (X(1)>0, 

X (N)  = 1.0) 

Y(I)  Velocity  ratio  change  at  X/RL  = X(I) 

COMMENTS  ON  USAGE 

1.  In  most  cases,  where  full  output  from  the  boundary  layer  program  is 
desired  only  for  the  final  pressure  distribution,  the  setting  of  1VF  <_  0 
will  give  a total  output  which  is  typically  one-third  of  the  corresponding 
output  for  IVF  >1.  In  the  case  when  IVF<0  and  IPF<0,  the  program  still 
prints  out  various  overall  boundary  layer  characteristics  at  the  tail  of 
the  body  as  well  as  the  error  in  C at  all  points  along  the  surface  of  the 
bodv-wake  displacement  model  for  intermediate  iterations. 

2.  For  input  values  of  NXT>NXMM  and  LGlb  = 0,  the  boundary  layer  is 
assumed  to  be  laminar  up  to  separation,  regardless  of  Reynolds  number  or 
pressure  distribution.  While  such  flows  arc  not  physically  meaningful,  they 
may  be  of  interest  in  certain  theoretical  studies  of  laminar  boundary  layers. 

(i 


3.  Since  the  boundary  layer  program  stops  at  separation,  the  program 
has  to  be  rerun  for  cases  where  the  separation  point  occurs  ahead  of  the 
input  value  of  XWAK,  which  typically  is  0.9S.  In  these  cases,  the  new  value 
of  XWAK  should  be  input  as  a value  just  less  than  TXl(I-2)/RL,  where  TX1(I)RL 
is  the  value  of  X/RL  where  separation  occurs. 

4.  As  pointed  out  above,  an  input  value  of  IPR0P>1  means  that  the 
velocity  distribution  is  not  calculated  by  the  program  but  is  instead  input 
into  the  program.  These  distributions  may  be  experimental  values  or  may  be 
the  results  of  a previous  calculation  using  this  program.  In  these  cases, 
the  program  skips  the  potential  flow  calculations  and  goes  directly  to  the 
boundary  layer  calculations. 

5.  When  IPR0P>1,  provision  is  made  in  the  program  for  reading  in 
changes  to  the  input  velocity  distribution.  These  changes  may  correspond, 
for  example,  to  modifications  of  the  input  bare  hull  velocity  distribution 
due  to  the  presence  of  an  operating  propeller  and/or  appendage.  Por  the  case 
where  no  changes  to  the  input  distribution  are  desired,  simply  set  the 
quantities  Y(I)  = 0 for  I = 1 to  N . 

6.  For  cases  when  a device  is  used  to  trip  turbulent  flow,  the  results 
of  Reference  8 indicate  that  the  tripping  device  has  measurable  parasitic  drag. 
The  net  effect  is  that  the  effective  location  of  transition  is  no  longer  at 
the  tripping  device  but  is  instead  moved  forward  to  a virtual  origin,  which 
depends  on  the  trip  location,  geometric  information  on  the  body,  and  computed 
laminar  boundary  layer  parameters  on  the  forebodv.  The  procedure  for  making 
this  calculation  is  given  in  References  1 and  8 and  hence  is  not  repeated  here. 
The  value  of  NXT  defined  above  should  be  set  equal  to  the  index  of  the  body 
station  at  the  virtual  origin. 

COMPUTER  PROGRAM  STORAGE  AND  PI ME  REQUIREMENTS 

On  the  COG  b700  computer  currently  in  use  at  the  Center,  the  program 
requires  a memory  of  approximately  145,000  octal  words  and  a period  of  70 
seconds  to  compile.  Program  execution  time  depends  on  a number  of  variables 
such  as  the  number  of  body  and  wake  points,  and  the  number  of  iterations 
required  to  arrive  at  the  final  results.  For  a typical  example  of  a bodv-wake 


displacement  model  described  by  170  points,  an  execution  time  of  approximately 
350  seconds  is  required  to  make  four  complete  calculations  of  the  pressure 
distribution,  the  boundary  layer  flow,  and  the  wake  flow  (ICUT  =3).  if  one 
uses  the  overnight  computer  priority  P2,  the  total  cost  for  this  run  is 
approximately  $60. 

PROGRAM  OUTPUT 

The  output  of  the  program  is  best  illustrated  by  means  of  a sample 
problem. 


SAMPLE  PROBLEM 

The  computer  program  will  be  used  to  calculate  the  pressure  and  boundary 

layer  on  an  axisvmmetric  body  designated  as  Model  4620-3.  The  body  is 

described  by  IS”  points  and  is  catalogued  on  permanent  file  CHHX46203TP60 , 

ID=CHHX.  The  body  has  a total  length  of  10.0  feet  and  a maximum  radius  of 

O^d’S  feet.  The  body  is  to  be  tested  in  the  wind  tunnel  at  an  air  speed  of 

217  feet  and  a corresponding  Reynolds  number  of  1.268  x 101 . A transition 

trip  wire  is  placed  on  the  forebody  at  X/RL  = 0.05.  With  the  procedure  outlined 

in  Reference  8,  the  virtual  origin  of  turbulence  is  calculated  to  be  at  X/RL 

= 0.015.  Transition  may  be  assumed  to  be  instantaneous.  The  iterations  for 

pressure  are  to  stop  when:  (a)  the  difference  in  the  pressure  coefficient 

0^  between  two  successive  iterations  is  less  than  0.01  at  all  points  along  the 

bodv-wake  displacement  surface,  or  (b)  the  number  of  iterations  for  pressure, 

after  the  initial  potential  flow  calculation,  is  equal  to  5.  (For  most  cases, 

the  maximum  difference  in  C is  less  than  0.02  after  three  iterations.)  Use 

P 

the  pressure  from  the  previous  iteration  to  compute  the  boundary  layer  and 
wake  flow.  Use  171  points  to  model  the  body-wake  surface.  Guess  that  the 
wake  velocity  reaches  0.99  of  free  stream  velocity  on  the  wake  surface  at 
X/RL  =1.15.  In  order  to  avoid  separation  in  the  stern  region,  linearly 
extrapolate  the  initial  potential-flow  velocity  distribution  to  the  tail  fo» 
X/RL>0.95.  Take  the  upstream  and  downstream  matching  points  for  the  5th-degree 
polynomial  to  be  at  X/RL  = 0.95  and  1.05,  respectively.  Use  the  short  printout 
options.  Velocity  profiles  are  desired  for  X/RL>0.90. 
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McCarthy,  J.H.,  J.L.  Rower,  and  T.T.  Huang,  "The  Roles  of  Transition, 
Laminar  Separation,  and  Turbulence  Stimulation  in  the  Analysis  of 
Axi symmetri c Body  Drag,"  Eleventh  Symposium  on  Naval  Hydrodynamics , 
London  (1976). 





SOLUTION 

The  data  cards  for  this  problem  are  listed  in  Table  1.  The  first 
card  attaches  the  permanent  file  describing  the  body  points.  The  other 
cards  follow  the  order  of  the  READ  statements  given  previously.  The 
symbol  b is  used  to  denote  a blank.  Also,  column  numbers  1,  11,  21,  31, 

41,  and  51  have  been  indicated  since  most  of  the  data  start  in  these 
columns . 

GENERAL  DESCRIPTION  OF  OUTPUT 

Since  the  short  printout  options  were  used  for  the  present  problem, 
the  program  prints  the  full  output  from  only  the  initial  pass  through  the 
potential  flow  program  and  the  final  pass  through  the  boundary  layer  and 
wake  programs.  Table  2 shows  sample  portions  of  each  section  of  the 
output.  The  portions  which  are  not  shown  are  simply  the  remainder  of  the 
station  data  for  each  section. 

The  potential  flow  program  first  prints  the  maximum  error  in  the 

source  density  for  each  Seidel  iteration  of  the  simultaneous  equations 

for  source  density.  Table  2 shows  that  a total  of  eight  iterations  are 

required  to  reach  the  convergence  criterion  for  the  maximum  allowable  error 

in  source  density,  whicii  is  1 x 10  * . The  potential  flow  program  then 

prints  various  geometric  and  flow  variables,  including  the  pressure 

coefficient  C , at  each  station.  After  the  initial  potential  flow  printout, 

P 

the  boundary  layer  program  prints  out  the  geometry  of  the  body  in  dimensionless 
and  dimensional  coordinates  as  well  as  other  geometry  and  input  variables. 

The  following  several  pages  of  Table  2 give  the  abbreviated  output  for 
intermediate  iterations.  The  boundary  layer  program  prints  out  the  difference 
in  C between  successive  iterations  at  each  point  along  the  body-wake 
displacement  surface.  Table  2 shows  that  this  maximum  difference  for 
iterations  1,  2,  and  3 is  respectively  0.235,  0.025,  and  0.013.  Hie  program 
also  prints  out  the  drag  of  the  body  calculated  by  three  separate  methods 
and  referenced  to  three  different  areas;  several  boundary- layer  variables 
at  the  tail  are  also  printed  out.  The  potential  flow  program  prints  out  only 
the  maximum  error  for  each  Seidel  iteration  for  source  density. 
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After  the  final  iteration  for  C , the  output  from  the  boundary- 
layer  program  is  greatly  expanded.  First,  the  pressure  coefficient 
for  the  initial  potential  flow  calculation  CPPF,  the  final  iterated 
viscous  pressure  coefficient  CP,  and  the  difference  CP-CPPF  are  printed 
out  at  each  station.  The  associated  body  pressure  drag  at  each  station, 
as  well  as  the  cumulative  total  pressure  drag  from  the  nose  to  the  station 
for  the  above  three  cases,  are  also  printed  out.  Detailed  velocity 
profiles  are  printed  next  over  the  rear  of  the  body.  These  velocity  profiles 
are  given  in  terms  of  components  normal  and  tangential  to  the  body  as  well 
as  in  the  axial  and  radial  directions.  Following  these,  a summary  of  local 
and  integrated  boundary  layer  variables  at  each  station  is  printed.  The 
program  concludes  by  printing  the  offsets  of  the  final  body-wake  dis- 
placement surface. 

DEFINITION  OF  OUTPUT  VARIABLES 

The  output  variables  are  defined  in  the  order  that  they  appear  in 
Table  2. 

Potential  Flow  Program 

\ Dimensional  axial  distance  in  feet  measured  from  the  nose  of  the 

body 

V Dimensional  radius  in  feet  of  body 

TL  Local  tangential  velocity/free  stream  velocity  UI 

■y 

CP  Pressure  coefficient  = 1 - (Tl) 

SIN  A Sin  a,  where  a is  defined  in  Figure  1 

COS  A Cos  a 

SIGMA  Source  density 

N Local  normal  velocity/UI 

PHI  Perturbation  potential  due  to  presence  of  body 

Boundary  Layer  Program  - Part  1 

TRTLAG  LG16  (see  DEFINITION  OF  INPUT  VARIABLES) 

TRINT  LG17  (see  DEFINITION  OF  INPUT  VARIABLES) 

TVC  = 1,  transverse  curvature  effects  are  taken  into  account 

SIIORTP  = 1,  velocity  profiles  are  printed  for  X/RL  > FLPRNT 


K Station  number 

X/C  Axial  distance/RL 

S/C  Arc  length/ RL 

Y/C  Radius/RL 

X Axial  distance  in  feet 

S Arc  length  in  feet 

Y Body  radius  in  feet 

RL  Body  length  in  feet 

RilOREF  p=  MURLF  * Rl/UI,  density  of  fluid  in  slugs/ft'5,  based  on 

input  Reynolds  number/ft  Rl  and  velocity  UI 

MUREF  y=  0.3834  x 10  b lb  sec/ft“,  dynamic  viscosity  of  air  at  60°F 

UI  Free  stream  velocity  in  ft/sec 

Rl  Reynolds  number  per  foot 

R1*RL  = R1*RL,  Reynolds  number 


Boundary  Layer  Program  - Part  2 


N 

Station  number 

X/C 

Axial  distance/RL 

S 

Arc  length  in  feet 

RO/C 

Body  radius/RL 

BETA 

B=  (2C/UE) (dUE/dt) 

term 

SQUIG, 


dimensionless  velocity-gradient 


CP  PF 


Pressure  coefficient  calculated  using  original  body,  i.e.,  C 
from  initial  potential  flow  calculation  ^ 


SQUIG 


R0  2 **2  4 

UE(— )“  ds,  transformed  arc  length  coordinate  in  lb^sec  /ft 


COS(ALPHA)  Cos  a,  see  Figure  1 for  definition  of  a 


SIN (ALPHA) 

CRINT 

UE 

CP 

mijf 


Sin  a 

= 2 AC  (R0/C)*tana/(R0MAX)‘ , ACp  = CP-CPPF 

Local  tangential  velocity  calculated  by  potential  flow  program 
using  latest  body-wake  displacement  model  in  ft/sec 
2 

= 1 - (UE/UI)  , pressure  coefficient  calculated  using  the  latest 
body-wake  displacement  model 

^ C ”1  Q 

= (1.3834  x 10"  lb  sec/ft*-,  dynamic  viscosity  of  air  at  60  F 
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CRLPF 


CRLVF 

CRLDF 

CRAPF 

CRAVF 

CRADF 

CRAFA, 

CRAWA, 

CRAV23 

Boundary 

I 

FTA 

F 

FP 

FPP 

Y 

EPS+ 

UPLUS 

WEL/UE 

X/C 

R/C 

R/RMAX 

R-RO/RM 


= 2 CP  PF  * DR/(ROMAX)  , CP  PF  = (CP  PF(N)  + CP  PF(N  - l))/2, 

DR  = R0(N)/C  - R0(N  - 1 ) /C , component  of  CP  PF  in  axial 
direction  referenced  to  frontal  area 

= 2 CP  * DR/(R0MAX)2,  CP  = (CP(N)  + CP(N  - l))/2,  component 
of  CP  in  axial  direction  referenced  to  frontal  area 

= CRLVF  - CRLPF 
N 

= E CRLPF(I),  cumulative  total  from  nose  to  station  N of  CRLPF 
1=1 

N 

= E CRLVF(I),  cumulative  total  from  nose  to  station  N of  CRLVF 
1 = 1 

N 

= E CRLDF(I),  cumulative  total  from  nose  to  station  N of  CRLDF 
1 = 1 

Drag  component  of  (CP  - CP  PF)  for  the  entire  body  referenced 
to  (frontal  area,  wetted  area,  (volume) 2/3) > note:  CRAFA  = 

CRADF (NXMM) 


Layer  Program  - Part  5 


Station  number 
rY 


n= 


[ (— )dy,  transfomed  y-coordinate , y 

Jo/2E  C 


>o/2E, 

normal  to  body  (see  Figure  1) 


distance  measured 


= - , 1 1>  = stream  function,  F = dimensionless  stream  function 

/n  (rl) 

= 8f/9n  = u/UE,  u = tangential  velocity  in  ft/sec 
= 3 2f/3n2  = 3(u/UE)/3n 

y,  distance  measured  normal  to  body  in  feet  (see  Figure  1) 
e+  =e/v  , e=  eddy  viscosity,  v=  u/p  = kinematic  viscosity 
u + = (u/UE)  / /c  J./2 , c}.  = local  skin-friction  coefficient 

= v/UE,  v = velocity  normal  to  body  in  ft/ sec 
Axial  distance  measured  from  nose/RL 
Radial  distance  measured  from  body  axis/RL 

Radial  distance  measured  from  body  axis/maximum  radius  of  bod) 

Radial  distance  measured  from  body  surface/RMAX , RO  = radius  of 
body  in  feet 
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ux 

* VUE'  “x  ■ 

UR 

■ UR/UE'  UR  ■ 

UTOT 

/ UX2  + UR2 

axial  component  of  velocity  in  ft/sec 
radial  component  of  velocity  in  ft/sec 


Boundary  Layer  Program  - Part  4 


N 

S 

THETA 

DELS 


CF 

CRL 

CD (CALC) 


Station  number 

Arc  length  measured  from  nose  in  feet 

9 = ) ttf-  O-nr)  dy,  momentum  thickness  in  feet 

J0  KU  Ufc  Ufc 

* P R u 

5 =1  — (1-— rr)  dy,  displacement  thickness  in  feet 

Jo  RO  Ut 

2 

C-  = t /(l/2pUE“),  r = shear  stress  at  the  wall,  local  skin 
r w w 

friction  coefficient  referenced  to  local  velocity  UE 

= CRLDF,  see  Boundary  Layer  Program  - Part  2 (BLPP2) 

= CRA(N)  «■  CFA,  CRA(h')  = CRADF  defined  in  BLPP2 , cumulative 

total  from  nose  to  station  N of  drag  coefficient  due  to  pressure 
and  friction,  referenced  to  frontal  area 


CD (GRAN) 


I MAX 

X/C 

RX 


RTHETA 
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4it(R0)  0 ,UeJ(H  + 2)  + 3^/8  nwAV  j.  r u j 

= — - — - — - (— ) , RMAX  = maximum  radius  of  body, 

*(RMAXr 

drag  coefficient  referenced  to  frontal  area  computed  by  the 
Granville  formula  from  nose  to  station  N 

Total  number  of  points  in  the  n-direction  required  to  reach 
the  outer  edge  of  the  boundary  layer  where  u/UE  = 1 

Axial  coordinate/RL 

= (UE)S/v,  v = kinematic  viscosity,  Reynolds  number  based  on 
arc  length  S and  local  velocity  UE 

= (UE)0/v,  Reynolds  number  based  on  momentum  thickness  9 and 

local  velocity  UE 

★ 

=6  /6,  shape  factor 


DEL (GRAN) 


CFA 


= + RO-  ♦ 2R0 1 cos  a ( effective  thickness  due  to  displace- 

cos  a 

ment  effect  of  boundary  layer  to  be  added  to  original  body  to 
obtain  body-wake  displacement  surface,  in  feet 


Cumulative  skin  friction  coefficient 
referenced  to  free  stream  velocity  U 


from  nose  to  station  N, 
and  frontal  area 


1 ’ 


CRA 


= CRADF  defined  in  BLPP2 


rnrcn  vmi  4tt(RO)0  ,UE , (H+5) / 2 , . 

CD(SQ-YN)  = — - — - — (yrr)  » drag  coefficient  referenced  to 

n(RMAXj  U1  5 

frontal  area  computed  by  the  Squire-Young  formula 

ETAINF  Value  of  n at  outer  edge  of  boundary  layer  where  u/UE  = 1 

CDT(CALC)  = CRAFA  + CRA,  CRAFA  is  defined  in  BLPP2,  sum  of  (total  drag 

coefficient  on  entire  body  due  to  pressure)  and  (cumulative 
skin  friction  coefficient  from  nose  to  station  N)  referenced 
to  UI  and  frontal  area 

WA.V23.FA  Subscripts  indicating  that  the  drag  coefficient  is  referenced 
to  (wetted  area,  (volume) 2/3,  frontal  area) 

N Station  number 

(X/C)OLD  Axial  distance  of  a point  on  the  original  body/RL 

RO  Radius  of  body  in  feet 

DELS  6*,  displacement  thickness  in  feet 

DELGR  = DEL(GRAN) , effective  thickness  in  feet  to  be  added  to 

original  body  to  obtain  body-wake  displacement  surface 

(RO  + DELG)/C  = (RO  + DELGR)/RL,  (body  radius  + effective  thickness) /RL , 
dimensionless  radius  of  body-wake  displacement  surface, 
assuming  the  effective  thickness  is  added  normal  to  the 
x-axis 

(X/C)NEW,  (RO  Value  of  (X/RL,  R/RL)  at  surface  of  body-wake  displacement  model, 

tr-  assuming  the  effective  thickness  is  added  normal  to  the  surface 
+ DG*COS)/C  - . . , , , 

of  the  original  body 

U TAIL,  X Value  of(UL/UI,  X/RL,  H)  at  tail  of  original  body 

TAIL,  II  TAIL 

">  rnr  f-a 

GAMA  TAIL  = (j^-)"  -~IL)  l7*11  TAIL  ♦ 2)  ♦ 3]/8,  initial  value  of 

dimensionless  momentum  area  for  the  start  of  the  wake  calculations, 
based  on  overall  drag  coefficient  CDC  FA  and  velocity  ratio  at 
the  tail  (derivation  is  given  in  Reference  1) 

Wake  Program 

X*  Axial  distance  measured  from  nose  in  feet 

NOM  AREA*  =\  R gp-  J>'  = (RO)  (0) , momentum  area  in  ft" 

FORM  FAC  H*  = DISr  AREA /MOM  AREA,  form  factor 

OD 

DISP  ARIA*  =\  R(l-mp)  dy  = (RO)  (6*),  displacement  area  in  ft" 

I n ut 
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m 


EFF  THICK 
11/ II  INF 
CF 


= / 2 ( D I SP  ARFA) , effective  radius  of  wake  in  feet 
= UE/Ul  in  wake 
= 1 - (U/UINI-T 

‘The  dimensional  definitions  are  given.  The  dimensionless  quantities  are 
divided  by  the  appropriate  power  of  body  length  RL . 

COMPARISON  OF  RESULTS  FOR  TWO  DIFFERENT  TRANSITION  LOCATIONS 

In  order  to  assess  the  importance  of  using  the  virtual  origin  as  the 
location  of  transition,  a second  computer  run  was  made  in  which  all  of 
the  input  parameters  were  kept  the  same  except  that  transition  was  assumed 
to  be  at  the  trip  wire  location,  X/RL  = 0.05.  Table  3 shows  the  results 
at  the  tail  of  the  body  for  the  two  assumed  locations  of  transition.  The 
differences  are,  on  the  whole,  not  large.  The  largest  differences, 
approximately  1 to  2 percent.,  occur  for  the  drag  coefficients  CIA  VVA, 

CRA  WA,  CDC  WA,  and  the  boundary  layer  displacement  thickness  DELS.  These 
are  due  to  the  parasitic  drag  and  boundary  layer  thickening  effect  of  the 
trip  wire.  On  the  other  hand,  the  shape  factor  H TAIL  and  the  velocity 
ratio  U TAIL  agree  to  nearly  four  significant  figures. 
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TABLE  1 

INPUT  DATA  FOR  SAMPLE  PROBLEM 
Column  Number 

1 11  21  31  41  51 

ATTACH , TAPE60, CHHX4 620 3TP60 , ID=CHHX. 

bbOl 37bbObbObO . 90  CALCULATIONS  FOR  MODEL  4620-3 

bbb400 

0 . 044  75  10.00  21 7 . 1268000. 

bl "lbbbSbbbO 

1.15  0.95  0.95  1.05  0.01 
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DOUGLAS  AIRCRAFT  COMPANY 
LONG  BEACH  DIVISION 

POTENTIAL  FLON  CALCULATIONS 


CP  ITERATION  1 

1 • J596AE-G I 

2 • 1 NN8  7E- 0 1 

i .J761JE-0Z 

<•  • AlObSE-O  3 

5 .13B15E-0J 

• 841E«00  .?78E»00  .687E«C3 

6 • 1 5 7N6E  -0  <• 

7 • 1087LE-05 

8 • 80  759E-06 

8 ITERATIONS  REQUIREO  FOR  CONVERGENCE 


N 

X MIDPOINT 

CP  e««o* 

1 

• C 0 25  0 0 

. 002320 

2 

.002500 

. 002431 

3 

.012503 

.002863 

4 

.017500 

.003427 

5 

.022000 

. 002507 

6 

. 026000 

.002155 

7 

. 030000 

.002033 

0 

. 033500 

• CO  1 990 

9 

. 030530 

. 001844 

10 

. 046000 

.001740 

11 

. 052500 

.001040 

12 

. 057500 

.001876 

13 

. 062500 

.001062 

14 

.067500 

.001007 

15 

. 072500 

. 0019*6 

16 

. 077500 

• 0 C 23  0 0 

1 7 

.082500 

.002094 

IS 

. 00  7500 

. 002173 

19 

• C925  0 C 

. 002218 

20 

. 097500 

. 002270 

21 

. 102500 

.002353 

22 

.137500 

. 002390 

23 

. 112500 

.002427 

24 

.117500 

. C02426 

26 

. 122500 

. C 0 23  06 

26 

. 127500 

. 002294 

27 

. 132500 

. 0021 39 

2 ft 

.137500 

. 001931 

29 

.142500 

. 001563 

30 

. 147500 

. C0U0  7 

31 

. 162500 

.000526 

J2 

. 1 5 7S0C 

.000109 

3 J 

• 162500 

. 00 062 3 

5** 

. 167500 

. 000064 

55 

. 1 72503 

.000073 

36 

. 1 7 7500 

.000791 

J 7 

. 102500 

.000703 

J0 

.107500 

. CQ06I2 

39 

. 192500 

. COO501 

HQ 

.197500 

• 004526 

HI 

.205000 

.000433 

Hi 

.215000 

• 000J79 

Hi 

•225000 

. OOOJ66 

HH 

.235000 

.000316 

H 5 

• 24.5000 

.000265 

46 

•255000 

• 0002bl 

Hf 

•265000 

.000256 

HQ 

.275000 

. 000241 

H 9 

•265000 

. 000229 

50 

•295000 

. 000230 

51 

•306000 

.000226 

52 

•315000 

. 000223 

53 

•325000 

.000221 

54. 

•335000 

• 000220 

55 

• 365000 

• 000219 

56 

.355000 

• 000219 

57 

•365000 

• 000220 

5« 

• 37  5 0 CO 

• 000221 

55 

. 3650C0 

•00CZ24 

60 

. 395000 

. 000227 

61 

• 4.05000 

. 000231 

62 

. 4.15000 

. 000236 

63 

. 4.25  000 

•000243 

64. 

. 4.J500C 

. COQi 51 

65 

• 4.4.50  00 

. 000261 

66 

.<•55000 

• 00C272 

67 

• 465000 

.000294 

6A 

.475000 

.000298 

69 

.465300 

.000314 

70 

. 495  0 C 0 

.000333 

71 

.505000 

• 0 C 0 3 55 

72 

• 5 1 5 0 0 C 

. CC0360 

73 

• c 25  0 0 0 

. 000439 

74. 

•535000 

. 00  3444 

75 

.545000 

. 000495 

76 

.555000 

. 000533 

77 

. 565  0 C C 

.000598 

78 

. 575  0 0 C 

. C0C637 

79 

. 5650C0 

. CCC586 

80 

.595000 

.000169 

61 

•605COC 

.CCC642 

62 

. 61 62  75 

. CO  1496 

63 

. 626362 

.C02077 

64 

.636893 

. 001674 

65 

. 645425 

.001793 

66 

. 654956 

.001542 

87 

.666394 

. CO  1434 

66 

.675926 

.031470 

69 

. 666456 

. 00133  6 

90 

.696695 

• 0C1268 

91 

. 706427 

• OC 1 359 

92 

.715959 

.G01J27 

93 

. 726491 

. OC  1439 

94. 

. 735022 

.001349 

95 

• 746460 

. 001453 

96 

. 755992 

. 001695 

97 

. 766523 

. 00  1647 

98 

.776961 

. 001604 

99 

. 766493 

. 0021  43 

100 

. 79t 024 

. 0C2226 

101 

. 605556 

.002596 

1C2 

. 6 1 1 0 87 

. CC2596 

1C  3 

. 826525 

. 002934 

104. 

• 8 3 to  J 5 7 

. CO  3624 

105 

. 645588 

. OG  3559 

(' 


10* 

• 69/02* 

. 00  3*8] 

10/ 

.666550 

. 00*921 

100 

.6/606* 

.00**5* 

10* 

.665621 

. 0061  )* 

110 

.6*5153 

.005660 

111 

. *065*0 

. 00  5*1* 

112 

.*16122 

. 002252 

11] 

• *2**66 

. 0061  38 

11* 

• 9J5060 

.003530 

US 

.**5000 

. 003026 

11* 

.*53500 

.00*580 

11/ 

.*5*500 

.015586 

116 

•*6*500 

. 02*25* 

11* 

. *66500 

. 03201  ) 

120 

.*/ 1000 

.038333 

121 

. *73000 

.0*5312 

122 

. 92*500 

. 0****2 

123 

.*25500 

.053306 

12s 

.*2/000 

.060*60 

125 

. *28500 

.066056 

12* 

.*80000 

. 02*651 

12/ 

.*82000 

.06*96* 

12* 

.*8*000 

.0*6322 

12* 

. *66 OCC 

. 106*2* 

130 

. *68000 

. 122*** 

131 

.**0000 

. 1 38280 

132 

.**2000 

.152001 

13] 

.**3500 

.122063 

13s 

.**5000 

. 1*1608 

135 

.**2000 

.2320*2 

136 

. 4*9300 

.23*526 

CP  ERROR 

.23*5/6 

iOtMhOS  iOdJDO 


OOUGLAS  AIRCRAFT  COMPANY 
LONG  BEACH  0IVISI08 


POTENTIAL  flop  calculations 


CP  ITERATION  2 

1 .3S9bbE-01 

2 .lbb87E-Ql 

] .37bllE-02 

b .B1050E-03 

5 .1380SE-03 

. 8b  If  * 0 0 . 277E  *00 

b .15720E-0N 

7 • 1091 bE -05 

8 • 8 0 b5  0E -0  b 

8 ITERATIONS  REQUIRED  FOR  CONVERGENCE 


H 

K MIDPOINT 

CP  ERROR 

1 

.002500 

. 000000 

2 

.007500 

.000031 

J 

.012500 

.000012 

4 

• 0 1 76  0C 

. 000037 

5 

. C 2 2 0 0 0 

. CCC00  3 

6 

. 02*000 

. CCOOOb 

7 

.030000 

.cocoon 

6 

. 03350(1 

.000001 

9 

.0  3S500 

.000001 

10 

. 056000 

.003001 

1 1 

.052500 

.000031 

1 2 

. 057500 

. 000302 

1 J 

. o 6 2 6 0 C 

.000031 

1 6 

• 0b 7560 

.003001 

15 

.072600 

.000032 

16 

.077500 

,00003b 

1 7 

.082500 

.00000b 

IS 

.067500 

. CC0009 

19 

.092500 

.000010 

2 0 

.697500 

. 0 t.  0 0 1 2 

2 1 

. 102500 

.00001b 

22 

.107500 

. 000021 

2 5 

. 1 12500 

. 00002b 

24 

. 1 1 7500 

. 0000  32 

2 5 

.122500 

. 0000  38 

c 6 

. 12  7500 

. OCCObb 

2 7 

. 1 32500 

. 000080 

28 

.137500 

. 000082 

29 

. 1425CC 

. C b 0 0 b 8 

3C 

. 167500 

. 0000  31 

31 

.152500 

.000008 

32 

. 15750  0 

. C0C083 

33 

. 1*2500 

.000030 

3 4 

• 16750C 

. CCC097 

35 

. 1 72500 

. 000077 

3 1 

. 1 77500 

. C000b9 

3 7 

. 1825CC 

. 000029 

Jrt 

. 1 6 7600 

.000017 

39 

. 19250C 

.000312 

. b 8 8E  ♦ 0 0 


r 


%Q 

. 191900 

.000007 

41 

. 209000 

.000002 

<.2 

• 2 1 S 0 00 

. ooooo i 

••a 

.229000 

.000001 

%% 

.219000 

. 000001 

Hi 

.299000 

. 000002 

Hb 

.299000 

.000003 

HT 

.269000 

.000003 

46 

•275000 

.000003 

*9 

.209000 

.00000% 

so 

.299000 

.000003 

SI 

. 309000 

.000001 

it 

. 319100 

. 000001 

53 

. 329000 

.000002 

5% 

. 339000 

.000003 

55 

.369000 

.000004 

56 

. 39900C 

. 000004 

91 

.369000 

.000004 

sa 

. 329000 

.000004 

ss 

. 389000 

.000009 

60 

. 395000 

.000009 

61 

.<•09000 

.000009 

62 

.415000 

.000005 

63 

.425000 

.000005 

6*» 

.<•39000 

.000006 

65 

.445000 

. cocao  6 

66 

.499000 

.000006 

67 

• 465000 

.000006 

66 

.419000 

.000006 

69 

.465000 

.000005 

70 

.495000 

.000036 

71 

.505000 

. 0 00004 

72 

. 515000 

.000002 

73 

. 525 DOC 

. DOQDDO 

74 

.939000 

. 000009 

75 

. 9490C0 

.000013 

76 

.999000 

. 0000  21 

77 

.565000 

• 000052 

76 

.929000 

. 000069 

79 

.565000 

.000107 

60 

. 596000 

.000062 

61 

.605000 

.000041 

82 

.616275 

.000166 

63 

.626362 

. 000225 

64 

.635693 

. 000163 

85 

• 645425 

. 0001  04 

86 

. 654956 

. U0Q062 

67 

.666394 

.000046 

66 

. 675926 

. 0000  46 

69 

. 665456 

. 0000  36 

90 

.696695 

.000035 

91 

. 706427 

• C0Q047 

92 

.715959 

. CG 0066 

93 

. 729491 

. C0C399 

94 

. 736022 

. oocosi 

95 

. 746460 

. COCO  70 

96 

. 756992 

. 0 0 CO  9 3 

97 

. 766523 

. 00Q0  86 

96 

. 776961 

. COCIO  7 

99 

. 766493 

.000146 

100 

. 796024 

.030184 

101 

.805566 

. 000204 

102 

.816067 

.000214 

103 

.626526 

. 0 00  10  3 

104 

. 636057 

• 0 0 C 4 3 1 

105 

. 645566 

.000433 

10b 

.057026 

.00060 3 

10  7 

. 066558 

• oooibb 

100 

. 076009 

• 0011 06 

109 

. 885621 

• 00 15b* 

no 

. 096153 

. 001950 

111 

. 906590 

. 003327 

lie 

.916122 

.005521 

113 

• 92b9b0 

. 007021 

lib 

.935000 

.00  779b 

ns 

. 9b  5 Q C 0 

. 007  225 

lib 

.953500 

.007002 

117 

. 959500 

.007633 

118 

. 96b  5 0 0 

. 006006 

119 

. 966S0C 

.006225 

120 

.971000 

. C 0 5bbb 

121 

.973000 

• 0 0 bb63 

122 

• 97b50G 

.003095 

123 

.976500 

.003321 

12*. 

. 977000 

.002223 

125 

. 97  050  C 

. OClbOb 

126 

.900000 

.000200 

12  1 

. 902000 

. 0C111 b 

120 

. 90*  0 00 

. CC2607 

129 

• 906000 

.00bl67 

130 

. 900000 

. 005702 

131 

. 990000 

.007637 

132 

. 992000 

.009110 

133 

. 993500 

. C 10665 

13b 

. 995 0C0 

. 0 1166b 

135 

.997000 

. 0 1 3355 

136 

. 9990  00 

. 015022 

137 

1 . 00  1000 

. C 15652 

130 

1 .003000 

.010236 

139 

1. 006000 

. 0 19766 

160 

1 . 00  7 0 CO 

. 021220 

161 

1 . G0900C 

. 02271 1 

192 

1 . 012000 

. 02  306.3 

1 b 3 

1 . 3170G0 

. C25390 

lbb 

1 . 026000 

. 0261  3 7 

lb5 

1 . 0 35000 

. 020167 

lb6 

l .0b5000 

. 011265 

1 b 7 

1 .05  7 000 

. CC1933 

lbO 

1 . 07200C 

.002666 

1 b 9 

1 . 09C0C0 

.002696 

150 

1 . 11GQ0C 

. 0 

151 

1 . 1 36  00  0 

. 00i7b7 

162 

1 . 1 75  00C 

.000960 

153 

t .22*  000 

.000610 

15b 

1 .275000 

. 000  3 75 

155 

1 . 325000 

. 000269 

156 

1 . J750C0 

. C 0 0 1 9 7 

157 

1 . bSCOOO 

.000121 

150 

1 . 5 5 C 0 0 0 

.000105 

169 

1 .675000 

. C0C077 

160 

1.07*300 

. 0bC0b7 

161 

2 . 200000 

. C0C026 

162 

2. 7CC00C 

. C jCQIQ 

1 6 J 

3.300000 

. C 0000  6 

16b 

b. 050000 

. 00300  1 

165 

5.25C0G0 

. G3CQ00 

166 

b.  750000 

. COCOGO 

167 

0. 75C0QO 

. 00000- 

160 

12.000000 

.000000 

169 

1 7.  COCOCO 

. cooooo 

170 

26.000000 

. CGOOOO 

MA  ■ CP  t WROP*  . j2539Q 


DOUGI.45  AIRCRAFT  COMPANY 
LONG  BEACH  OlVISIOR 


POTENTIAL  FLON  CALCULATIONS 


CP  ITERATION  J 

1 • J596NE-0 1 

2 .H.i.aFE-01 

3 .37611E-07 

<•  .81053E-03 

5 .1JA89E-03 

• SNLf»0:  .Z79f*0C  .fcAAf*:; 

6 .1S726E-0A 

7 .10906E-05 

» .8056AE-06 

# ITERATIONS  RE  QU  I RE  0 FOR  CONVERGENCE 


H 

X NIOPC1NT 

CP  ERROR 

1 

.002500 

.000000 

2 

•007500 

.000001 

3 

.012500 

.003001 

9 

. 017500 

.000001 

5 

.022000 

.000001 

6 

.026000 

.000001 

7 

. :3lOQC 

.000001 

a 

.033500 

. 00000  1 

s 

.036500 

.000001 

10 

• 096000 

.000331 

11 

.052500 

• C0G001 

1 2 

• G 5 75  C G 

.003001 

13 

. 062500 

. 000031 

1 9 

.067500 

. 0 0000  1 

15 

. 072500 

. 0C0001 

lb 

• C 77500 

.003001 

17 

.082500 

.000301 

18 

.087500 

.000001 

19 

.092500 

. C 000C  1 

20 

• 0975C0 

.000001 

21 

.102500 

• 0 w JO  3 2 

22 

. 10  7500 

. 300002 

23 

.112500 

. CC0002 

2* 

.117500 

. C 000  3 2 

25 

. 1225C0 

.000033 

2b 

.127500 

. CQQQQ 3 

27 

.132500 

• C C C 0 0 9 

28 

.137500 

. CCC005 

2 9 

.192500 

. CCGOOb 

30 

.197500 

.000305 

31 

. 152500 

.000002 

32 

.157500 

. OCOOO  3 

33 

• 16  250  3 

.000307 

3*. 

. 1675CG 

. C0C30  7 

35 

. 172500 

. 0 0 G 0 0 9 

3b 

.177500 

. OOCOOO 

37 

.182500 

. 00000  3 

38 

. 1 8 7500 

. OGOOO  3 

39 

.192500 

. COCO  33 

33 


%0 
41 
HZ 
Hi 
HH 
H 5 
46 
*.7 
48 
4S 
>0 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
b 7 
68 

69 

70 

71 

7 2 

73 

74 

75 

76 

77 

78 

79 

8 C 
81 
82 

83 

84 

85 
8b 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
lul 
K2 
i:  3 

1G4 

105 


s 


. I 97»oo 

• 205000 
.715000 

• 72  5000 
.235000 

• 7*5000 
.755000 
.765000 
•775000 
.705000 
.795000 
. 305000 
. 315000 
. 375000 
. 335000 
. 355000 
. 355000 
.365000 
. J750CJ 
•305000 
.395000 
.405000 
.415000 
.475000 
. 435000 
.445000 
.455000 
.465000 
.475000 
.485000 
. 495000 
. 505000 
. 5 1 5 0 0 4 
.575000 
.535000 
.545  000 
.555000 
.5e50C0 
.575000 
.585000 
.595000 
.605000 
. 6167  75 
.676367 
. 635893 
.645475 
.654956 
.666394 
.675976 
.685458 
. 696895 
. 70  t 47  7 
. 715959 
.775491 
. 73507? 

• 746460 
. 75599? 
. 765573 
. 776961 
. 786493 
. 796074 
. 805556 
. 815087 
.876575 
. 836057 
. 84558# 


. 000003 
.000003 
. 000007 
.000007 
.000007 
. 00000? 
.00000? 
.000007 
.000007 
.000007 
.00000? 
.00000? 
.000007 

. tcooo? 
.00000? 
.000003 
.000003 
. 00000  3 
.000003 
.000003 
.000003 
.000003 
. 00000  3 
.000003 
.000003 
.000003 
. C0C004 
. 0 3 4 0 0 4 
.000004 
.000004 
.000005 
. -OOOJ5 
• CO  0 0 3 5 
.000006 
.000006 
.000006 
.000004 
. CCOOOl 
.000010 
.000014 
. C03008 
.000010 
.000077 
. 0000  37 
.000073 
. CCOOC  9 
.000003 
. 0 0 000  3 
.000004 
. 000004 
.000005 
.000007 
.000008 
.000007 
.000005 
. 000006 
. 00000  7 
.000037 
•000000 
.000007 
.000030 
.000011 
■ 0000  7 3 
. COCO  33 
.400051 
.000088 
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TABLE  3 


COMPARISON  OF  RESULTS  AT  TAIL  FOR  TWO 
LOCATIONS  OF  TRANSITION 


Virtual  Origin  Trip  Location 

X/RL  = 0.015  X/RL  = 0.05 


CFA  WA 

0.002928 

0.002896 

CRA  WA 

0.000102 

0.000104 

CDC  WA 

0.003030 

0.003000 

DELS 

1.5914 

1.5734 

H TAIL 

1 .2765 

1 .2768 

U TAIL 

0.9490 

0.9489 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

111  DTNSRDC  REPORTS  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE.  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(2)  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE.  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(31  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS  NOT  FOR  GENERAL  DISTRIBUTION 


